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ABSTRACT 

Due to the large kick velocities with which neutron stars are born in 
supernovae explosions, their spatial distribution is more extended than that of 
their progenitor stars. The large scale height of the neutron stars above the disk 
plane makes them potential candidates for microlensing of stars in the Large 
Magellanic Cloud. Adopting for the distribution of kicks the measured velocities 
ON ; of young pulsars, we obtain a microlensing optical depth of r ~ 2A 10 X 10~ 8 

(where N w is the total number of neutron stars born in the disk in units of 10 10 ). 
The event duration distribution has the interesting property of being peaked 
at T ~ 60 - 80 d, but for the rates to be relevant for the present microlensing 
searches would require N w ^ 1, a value larger than the usually adopted ones 
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(JVio~0.1-0.2). 

CO 

The ongoing searches of baryonic dark matter in the Galaxy by means of microlensing 
of stars in the Large Magellanic Cloud (LMC) have produced very surprising results 
recently. Indeed, the first five events obtained by the EROS (Aubourg et al. 1993) and 
MACHO (Alcock et al. 1995) collaborations, with durations T ~ 10 - 29 d, suggested 
the presence of compact objects with masses ~ O.1M , i.e. near the brown dwarf range. 
However, the new analysis of the first two years of MACHO observations (Alcock et al. 
1996) has now returned events with longer durations (besides eliminating their two shortest 
duration events). These new results are indicative of larger lens masses, in the range 
0.1 - 1M Q , being responsible for the events. In addition, the large event rates observed 
require a considerable amount of compact objects to be present in the Galaxy. 

A first difficulty associated with the new range of inferred masses is that the lensing 
population can no longer be completely ascribed to the unseen brown dwarf counterpart 
of one of the Galactic stellar components, such as the thick disk (Gould, Miralda-Escude 
and Bahcall 1994) or the spheroid (Giudice, Mollerach and Roulet 1994). Although a 
population of white dwarf remnants in these components could help to fit the observed 
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MACHO event durations, it would require an initial mass function quite contrived in order 
to have ^ 90% of the mass of these components ending up in white dwarfs, and even in 
this case the resulting rates would still fall short with respect to the observed ones. On the 
other hand, the proposed interpretation of the lenses being an old white dwarf population 
in the dark halo makes the problem of the absence of any luminous stellar counterpart more 
severe than for the previous scenarios. 

In this letter, we want to study the relevance for the microlensing searches of the 
other natural candidate for dark lensing objects in the Galaxy: the neutron stars (NS). 
An important point to take into account in the study of the population of neutron stars in 
the Galaxy is the fact that the violent birth of these objects in supernova (SN) explosions 
makes them acquire large 'kick' velocities, as evidenced by the observed velocities of young 
pulsars (fastly rotating magnetized neutron stars). The pulsar velocities have been studied 
by Lyne and Lorimer (1994), who obtained a fit to the observed distribution of transverse 
speeds v t given by 

£0.13 

where t = v t /(330 km/s). Under the assumption that the kick distribution is isotropic, 
the two dimensional distribution is just a convolution of the underlying three dimensional 
distribution /^(v), 
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f2D(Vt)=V t l dv fMv) ■ (2) 
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Hence, / 3 £> can be reconstructed by solving an Abel integral equation, to obtainF] 
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/sdO) = -- / dv t d j 2D r^^ =- (3) 

7T Jv dV t , /„2 _ v 2 



The average velocity resulting from this distribution is (v) ~ 390 km/s. This implies that 
the NS born from a given Galactic population (e.g. the disk or the bulge) will end up with 
a completely different, and much more extended, spatial distribution. The physical origin 
of these kick velocities is still unknown, and proposed mechanisms include the unbinding 
of a close binary system (Dewey and Cordes 1987, Bailes 1989), an asymmetric explosion 
(Burrows et al. 1995) or the interaction of the emitted neutrinos with the NS magnetic field 
(Kusenko and Segre 1996). 

The spatial distribution of the NS born from massive disk stars during the whole 
lifetime of the Galaxy has been considerably studied as a possible Galactic source of gamma 



The result obtained here is equivalent to the one of Terasawa and Hattori (1995). 
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ray bursts (Shklovskii and Mitrofanov 1985, Li and Dermer 1992, Podsiadlowski, Rees and 
Ruderman 1995). The neutron stars density can be obtained from Monte Carlo simulations 
of the initial distribution of neutron star locations and kick velocities at birth, and then 
following their orbits in the Galactic gravitational potential up to the present (Paczyhski 
1990, Hartmann, Epstein and Woosley 1990, Terasawa and Hattori 1995). In particular, 
in order to fit the isotropy of the BATSE observations (Meegan et al. 1992), the gamma 
ray bursts in these Galactic scenarios must be produced by the very distant tail of the 
spatial NS distribution, i.e. at Galactocentric distances r ^ 100 kpc. As we will show, 
a large fraction of the neutron stars born in the disk (~ 30%) contributes to the density 
at radii 8 kpc < r < 40 kpc (having a considerable height above the disk plane). Hence, 
independently of their possible relevance to account for the gamma ray bursts, the NS can 
effectively act as lensing objects in the microlensing searches in the direction of the LMC. 

To compute the NS distribution we follow closely the work of Paczyhski (1990), 
but with the currently accepted velocity distribution, eqs. ([HI). We take for the radial 
distribution of the initial birth positions 

dP(R) oc exp(-R/h R )RdR, (4) 

taking for the disk scale length h R = 3.5 kpc. For simplicity, we take z= for the initial 
NS vertical positions, since anyhow very massive disk stars form with a small scale height 
(< 100 pc), so that this assumption is not essential for the microlensing results. With these 
distributions, we produced 1.6 x 10 5 initial values of positions and velocities and followed 
the orbits by integrating the equations of motion (using cylindrical coordinates) 

+ |, (5) 
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in the Galactic gravitational potential In addition, we have the equations stating the 
conservation of the angular momentum and the energy 
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J? = Rvrh = const 



E = V R + V z + V 4> + $ ( jR ' Z ) = COIlst - ( 7 ) 

In particular, we used the energy conservation to check that the final accuracy for this 
quantity after the numerical integration was better that 10~ 6 . 



We model the Galaxy with three dynamical components, a disk, a bulge and a dark 
halo, with $ = $^ + ${, + For the disk and bulge components we used Miyamoto and 
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Nagai (1975) potentials given by 
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(i = d,b). 
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For the halo we took a logarithmic potential (Binney and Tremaine 1987) 
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+ const, 
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allowing in principle for a non-spherical distribution (q^ ^ 1) in order to explore the 
effects of a flattened halo. The corresponding mass densities can be obtained from the 
Poisson equation V 2 $ = 4irGp. The parameters in eqs. (H) are chosen so as to produce an 
acceptable rotation curve, and are M d = 8.07 x 1O 1O M , = 3.7 kpc, q = 0.2 kpc for the 
disk and M b = 1.12 x 1O 1O M , a b = 0, c b = 0.227 kpc for the bulge, with v = 200 km/s and 
R c = 10 kpc for the halo potential. 

The resulting local column density, for a solar Galactocentric distance of 8.5 kpc, 
is So(l z l < 11 kpc) = 79Af /pc 2 (for q^ = 1), which is consistent with the dynamical 
estimates from vertical motion of stars (Kuijken and Gilmore 1991), and increases to 
~ 85M /pc 2 for q<p = 0.8 (which corresponds to and axis ratio in the density distribution 
q p ~ 0.4 (Binney and Tremaine 1987)). The rotation curve velocity, t> 2 = Rd^/dR | 2= o, is 
clearly insensitive to q^. 

Regarding the total number of neutron stars produced during the whole Galaxy 
lifetime, this is a quantity which is not easy to estimate, since it depends crucially on the 
initial mass function of heavy stars (m > 8M ). This is quite uncertain, mainly because 
of the unknown evolution of the star formation rate from the birth of the Galaxy up 
to our days. If the present accepted rate of Type II and lb (core collapse) supernovae 
of ~ 0.02 — 0.03 yr _1 , which is consistent with the estimated pulsar birth rate (Lyne, 
Manchester and Taylor 1985), were to be representative of the average SN rate in the past, 
one would estimate a total number of neutron stars produced in the disk N ~ 3 x 10 8 . 
However, the star formation rate was certainly higher at earlier times, mainly due to the 
larger fraction of gas present in the past and to the possibility of an initial burst of star 
formation. This can enhance a lot the fraction contributed by the heavy stars to the initial 
mass function (Larson 1986). For instance, a recent modeling of the SN production in the 
Galaxy (Timmes, Woosley and Weaver 1995, 1996) concluded that a more reasonable value 
for the total number of neutron stars is N^s — 1-9 x 10 9 , but the uncertainty in this number 
is undoubtedly large. In view of this, we will not adopt a definite value for N NS , giving the 
results in terms of N w = iV^rs/lO 10 , hoping that in the future this number will be known 
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more precisely. Since in the models leading to N 10 ^0.1 most of the SN explosions took 
place in the early life of the Galaxy, we just assumed for the Monte Carlo simulations that 
the NS were all born 10 10 years ago. The results are in any case nearly unchanged for a 
uniform distribution of birth times. 

Figure 1 shows the number density contours of the resulting NS distribution. This 
distribution is more extended than the one obtained by Paczyhski (1990) and Hartmann, 
Epstein and Woosley (1990), due to the smaller velocities adopted at the time by those 
authors. We have indicated also in the plot the coordinates of the line of sight to the LMC 
(which position at (R,z)lmc — (42, —26.5) kpc is also labeled), assuming a distance to it 
of D os = 50 kpc. Looking at this curve, it is useful to recall that using the parameter 
x = D i/D os to describe the normalized distance along the line of sight to an hypothetical 
lensing object, at distance D Q i from the Sun, one has x =z/zlmc- 

From the resulting NS density we can now compute the microlensing optical depth (for 
reviews on microlensing see Roulet and Mollerach (1996) and Paczyhski (1996)) 

r = / dx— : 
Jo ax 

with 

dr 4ttG ^ 9 

te = -#- I % x 0--z)p(z), (n) 

where p is the mass density of the lenses, for which we just assume a common NS mass of 
1.4M . The resulting value is r = 1.9N 10 x 10~ 8 . Also, defining 

^ = U ldx ^ D "- (12) 

we find an average distance to the lenses of (D Q i) ~ 8.8 kpc, which shows that the lenses 
are indeed quite far away (we recall that (D ol ) ~ 1.1,3.6 and 14 kpc for objects in a thin 
disk, thick disk and dark halo respectively). 

Another important quantity to describe the microlensing events is the distribution of 
event durations. To obtain it one needs the distribution of the velocities with which the 
lensing objects cross the line of sight to the LMC. To this end, we obtained the velocities of 
the NS which were near to the line of sight, computed the relative velocities with respect to 
this line, and then took the projections of these relative velocities orthogonal to the same 
line. If v D and v s are the velocities of the Sun and the LMC respectively, the motion of the 
line of sight (and hence of the so-called microlensing tube) is just v t = (1 — x)v Q + xv s . In 
a coordinate system where x points toward the Galactic center and y along the direction 
of increasing longitudes, we adopted, following Griest (1991), v D = (9,231, 16) km/s and 
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v s = (53, —160, 162) km/s. The quantity of interest is 

v± = |v r - (v r • L)L\, (13) 

where v r = v/ — v t is the relative velocity, with v/ being the lens velocity which we obtained 
with the Monte Carlo. L = (cos b cos £, cos b sin £, sin b) is the versor in the direction of the 
LMC, which has Galactic coordinates (b,£) LM c — (—33°, 281°). 

The duration of the events is given by T = Re/v±, with the Einstein radius being 
Re = 2y ' GmD os x{\ — x)/c. The inclusion of the tube velocity is relevant in this scenario, 
contrary to the case of lenses in a halo, mainly due to the importance of the observer's 
motion and the 'memory' that the NS have of their progenitor's motion of rotation. 

We can compute now the 'theoretical' rate of events (assuming unit efficiency) from 
where 

dx ndx{T)(xY 1 ' 

with (T)(x) = R E (x)(vl 1 )(x). The resulting value is (T) th = 0.4iVio events/(10 7 stars yr). 
We can equally compute the 'theoretical' average event duration from the relation 



(T) = ^Tr^' (16) 

7T (1 }th 



resulting in (T) ~ 115 d. 



The event duration distribution can be obtained from the underlying distribution of 
velocities v± as a function of x. We have computed it by dividing the event durations 
in bins of AT = 10 d, i.e. considering the intervals [Tj,Tj + AT] with Tj = AT x i, 
and obtaining the fraction fi(x) of the events with durations in the specified intervals to 
estimate the differential distribution df{x,T)/dT of this fraction. The differential rate can 
then be computed using 

dT 2 /-i dr\ df(x,T) 



(T) = - / dx— - J \' ' . (17) 
dT n Jo dxT dT 

The resulting distribution is shown in Figure 2. A distinctive feature of it is that it peaks 

at T ~ 70 d, a fact that is interesting in view of the long duration of the events recently 

obtained by the MACHO collaboration (Alcock et al. 1996). However, one has also to 

recall that for their accumulated statistics of ~ 1.8 x 10 7 stars yr, and adopting an average 

efficiency in this range of event durations e ~ 0.3, one would expect a number of events 

N ev ~ 0.3iVio arising from the NS. 
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We have also explored the sensitivity of the results to variations in the assumptions 
made. Considering a flattened halo changes little the previous predictions. This is because 
the effects of flattening only affect sizeably the distribution of NS along the line of sight 
for |z| > 10 kpc, and the contribution to the optical depth of lenses in this region is not 
dominant in any case. 

Variations in the radial distribution of initial positions can have some effects on the 
predictions. For instance, adopting a Gaussian distribution centered at 5 kpc with an 
extension of 2 kpc, as done by Hartmann et al. (1990), increases the optical depth by 
~ 20%. Considering NS produced in the inner region of the Galaxy, as e.g. resulting from 
SN explosions of stars in the Galactic bulge, leads to a smaller optical depth 'per star 
borned'^, with events having somewhat larger durations. 

There is little sensitivity of the predictions to the details of the assumed velocity 
distribution, and for instance halving all the initial velocities leaves almost unchanged both 
r and (T). Most of the contribution to r comes from NS with intermediate velocities 100 
km/s & v £s 500 km/s. For example, if all stars were given initial kicks of 50, 100, 300 or 
600 km/s, we would obtain an optical depth r = 0.5, 1.6, 3.0, 0.5 x 10~ 8 iVio respectively. 

Finally, we would like to note that another logical possibility is that a fraction of the 
SN explosions leads to black hole (BH) formation rather than to NS. For instance, Timmes 
et al. (1996) have estimated that if all stars more massive than 19M gave rise to BH 
rather than to NS, ~ 70% of the remnants would be BH (Nbh — 1-4 x 10 9 according to 
their estimates). Two important differences would result in this case. The kick velocities 
imparted to the BH in the explosions, if any, would certainly be smaller, leading to a BH 
distribution less extended than the one obtained for the NS, and this would tend to reduce 
the optical depth toward the LMC. Second, the larger average masses of the produced BH 
would tend to enhance the optical depth proportionally (r oc p oc (ttibh))- If a significant 
fraction of the BH acquire velocities larger than ~ 100 km/s, the first effect will not be 
large, as previously discussed, so that the second one would dominate, allowing for an 
increased value of r. The larger masses would also imply longer event durations (T oc \fm). 

In conclusion, we have considered the contribution that NS may provide to the 
microlensing of LMC stars. Leaving aside the possibility of a pregalactic (population III) 
NS component (Eichler and Silk 1992, Brainerd 1992), the largest NS population would be 
the one arising from SN explosions in the Galactic disk, since this is the most massive stellar 
Galactic component. The resulting optical depth for this population is r ~ 2 x 10 _8 iVio, 
and the event durations are peaked at T ~ 70 d. This range of durations is interesting in 



2 For instance, r ~ 0.8 x 10 8 iVio for NS born in the inner 1 kpc of the Galaxy. 
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view of the long event durations observed by the MACHO collaboration, with the longest 
two having T = 65 and 50 dQ (a larger one, with T = 71 d and produced by a binary lens[], 
most probably corresponds to a lens in the LMC, due to its measured proper motion, while 
the events of shorter durations may (at least partially) be ascribed to a faint component of 
some of the Galactic or LMC stellar populations (see de Riijula et al. 1995)). However, 
for the reference value N w ~ 0.2, the expected event rate from the NS falls below the 
observed one by an order of magnitude (although the statistics is certainly not large). Only 
for a larger NS density, iVio ^ 1, the rates from NS would be significant for the present 
microlensing searches. Whether such large values are acceptable and consistent with the 
Galactic chemical evolution remains to be seen. At any rate, the existence of this lensing 
population certainly made it deserve a quantitative study. 



We thank A. Lanza for useful discussions. 
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Fig. 1. — Present neutron star density contours vs. cylindrical galactic coordinates. The 
different contours correspond to N/N w =3.E-3, l.E-3, 3.E-4, l.E-4, 3.E-5, l.E-5 (stars/pc 3 ). 
The location of the LMC is indicated as well as the coordinates (solid line) of the line of 
sight to that galaxy. 



Fig. 2. — Differential event rate distribution vs. event duration, taking as normalization 
iVio = 1. 
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